ANALYSIS OF A MATHEMATICAL MODEL FOR THE 
GROWTH OF CANCER CELLS 



MARTIN KOHLMANN 

Abstract. In this paper, a two-dimensional model for the growth of multi- 
layer tumors is presented. The model consists of a free boundary problem for 
the tumor cell membrane and the tumor is supposed to grow or shrink due 
to cell proliferation or cell dead. The growth process is caused by a diffusing 
nutrient concentration a and is controlled by an internal cell pressure p. We 
assume that the tumor occupies a strip-like domain with a fixed boundary at 
y = and a free boundary y = p(x), where p is a 2-Tr-periodic function. First, 
we prove the existence of solutions {a, p, p) on a scale of small Holder spaces 
and show that our model allows for flat stationary solutions. As a main result 
we establish that these equilibrium points are locally asymptotically stable 
under small perturbations. 
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1. Introduction 

The mathematical modeling of cancer growth is a challenging area of research 
in the applied sciences nowadays. The complex growth process of a tumor cell 
can be captured within different mathematical models, e.g., models consisting of a 
system of coupled partial differential equations which arise from reaction-diffusion 
equations and a mass conservation law 0, 0, El[ 28 , 2^] . In this context, tumor 



growth is often considered as a free boundary problem 0, Q, 13-Eol 12-15, Ec 
Sometimes, it is also justified to treat cancer growth as an incompressible flow 
in a porous medium so that cells move in accordance with Darcy's law, see, e.g., 
[si El, ■ We refer the reader to the review papers [1, 0, which present a 
variety of other tumor growth models. 



2000 Mathematics Subject Classification. 35R35, 35Q92, 35B35. 

Key words and phrases. Free boundary problem, tumor growth, well-posedness, flat stationary 
solutions, stability. 

1 



Analysis of a mathematical model for the growth of cancer cells 



2 



Over the years, the following general non-dimensionalized moving boundary 
value problem 



(1) 



Act = f(&) in 

Ap = -g{a) infi(t), 

a = a on dft(t), 

p = 7>t on d£l(t), 

V = -d v p ondfl(t), 

has turned out to be very appropriate to describe the growth of tumor cells; see, 
e.g., [l9| . Here a = a(t, x) and p = p(t,x) are unknown functions defined on 
the time-space manifold U t>0 {i} x tt(t) with an a priori unknown time-dependent 
domain Q(t) C R". The boundary dft(t) has to be determined together with the 
functions a and p which model the nutrient concentration and the internal pressure 
for the tumor cell described by D,(t). Furthermore, er > is a constant and k and d v 
denote the mean curvature and the normal derivative (with respect to the outward 
normal) for the boundary 9f2. Finally, 7 > is the surface tension coefficient 
and V denotes the normal velocity of the free boundary. In (HJ, the tumor f2(i) 
receives a constant supply of nutrient on its boundary and the pressure on d£l(t) 
is proportional to the curvature of the free surface dfl(t), as proposed in fjjj ]. The 
evolution equation for the moving boundary comes from an application of Darcy's 
law from fluid mechanics. We also impose the initial condition O(0) = Q,q, where 
f2o C 1™ is a given bounded domain in M™ with a sufficiently smooth boundary. 
Typical choices for / and g are 

/(er) = Act, g(a) = (t(a - a) 

where A, /i,<r > are constants, 0], or 

Aa mi 



/(*) = P nn , _ mi +H°), 



Act"* ( 5a— 



with positive constants A, B, ft, S, mi, m 2 , cr c , and a non-negative increasing 



function h, [281 ] . In the authors assume that / and g are general functions 
satisfying 

. f,g eC°°[0,oo), 

• f'{a) > for cr > and /(0) = 0, 

• g'(a) > for a > and g(a) = for some a > 0, 

• cr < a. 

In this paper, we will suppose that 

f{°) = V, ffW = m(ct - cr), 

with positive parameters fi, a. Since Ap is minus the divergence of the cell velocity 
field, the meaning of Ap = — /i(cr — cr) is that tumor volume is produced, if a is 
above the proliferation threshold a, and that the tumor volume decreases in the 
opposite case. 
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The tumor growth model which we study in this paper has the following form: 
The tumor is assumed to occupy a two-dimensional region of the form 

n,(t) ~{(x,y)eR 2 ; 0<y<p(t,x)} 

where t is the time variable and p(t, x) is an unknown positive 27r-periodic function. 
The upper boundary of the tumor is denoted as 

r p (t) := {{x,y)eR 2 - y = p(t,x)}, 

its lower boundary is To = {(x, y) S R 2 ; y — 0}. A similar situation is discussed in 
[ll ] where the authors explain that the strip-shaped model refers to the growth of 
multi-layer tumors, a kind of in vitro tumors cultivated in laboratory by using the 
recently developed tissue culture technique, 23l 124.127 1. While there are only a few 



works dealing with strip-shaped domains, a variety of papers considering radially 
symmetric models for tumor growth have been published, cf., e.g., the seminal 
paper 



!J 



T p (t) 



Figure 1: A free boundary problem modeling multi-layer tumors. 



The following conditions on the tumor growth process are imposed: We assume 
that the tumor is constantly supplied with nutrient on dfl p (t); precisely, the nutrient 
concentration is 02 > on T p (t) and u\ > on IV As in ([lj, we assume that the 
pressure on T p is given by 7«r p ■ Since tumor cells should only grow in the positive 
y direction, we assume that p y = on IV Thus we are led to study the following 
system of equations: 



(2) 



Act 


= a 


in Q p (t), 


t > 


Ap 


= -/i(o--cr) 


in Q p (t), 


t > 


Pi 


dp 

~ dv 


on T p {t), 


t > 


(7 


= &2 


on T p (t), 


t > 


a 


= (Tl 


on r , 


t > 


P 


= 7^r p 


on T p (t), 


t > 


Py 


= 


on r , 


t > 


P 


= Po 


for t = 0, 





where v = (— p x , 1) denotes the outward normal on T p (t) with respect to fl p (t) 
and po is a given periodic function. The third equation in ([2]) can be derived as 
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follows: we assume that the normal velocity V of the boundary T p is equal to the 
cell movement velocity in the direction of the outward unit normal 

1 



Vq = 



pi 



cf. [EH, ESS El- Then Darcy's law and the relation V = p t (l + p 2 x )- 1/2 (cf. Q) 
imply that V = — vq ■ Vp and that the motion of the free surface T p is modeled by 
the third equation of ©. 

In the model presented in [12j , the lower boundary Tq is supposed to be imper- 
meable for glucose and oxygen. The authors also comment briefly on a variant of 
their model obtained by exchanging the boundary conditions for a on T p and Tq. 
The novel aspect of the problem presented in the paper at hand is that we allow 
for supply of nutrient on both boundary components. 

Our paper can be outlined as follows: In Section [2j we recall some elementary 
facts, definitions and notation from [121 ] which will be important for our approach 



to the system ([2]). In Section [3l we prove that the system ((2|) is well-posed on a 
scale of small Holder spaces and we compute its flat stationary solutions. Finally, 
in Section [H we linearize the system @ at such an equilibrium and prove that it 
is asymptotically stable under small perturbations. 



Acknowledgement. The author thanks the anonymous referees for helpful sug- 
gestions that improved the preliminary version of the paper. 

2. Preliminaries 

Let § = M/27rZ. Function spaces of 27r-periodic functions will be identified with 
the corresponding spaces on S. In the following C+([0,T) x S) stands for the cone 
of all positive functions in C([0,T) x S), for any T > 0. We will make use of the 
little Holder spaces h k+a (§) which are defined as the closure of C°° (§) in the usual 
Holder space C k+a (§), for k G N and a G (0, 1). Similarly, for any open set C R 2 , 
we define h k+a (0) as the closure of C°° (H) in C k+a (0) . The little Holder spaces are 
Banach algebras under pointwise multiplication and the embedding /i r (S) h s (E>), 
for r > s, is compact. Wc shall denote the cone of positive functions in h k+a (E>) by 

/4 +a (§). 

Let 

D p . T = {(t, x, y); t G [0, T), x G S, < y < p(t, x)} . 
We call a triple (a, p, p) a solution to the problem ^ if 
(a,p) G C(D PtT ) x C(D P , T ), 

a(t, ■) G /i 4+a (p)), p{t, •) G h^iQplt)), t G (0,T), 
P C([(D) T), ft,^ +a (§)) n C 1 ((0, T), /i 1+Q (§)), 
[o~,Pip) satisfies © pointwise on D p t- 

A solution (a,p, p) is called stationary, if (cr,p, p) = {a{x, y),p(x, y), p(x)). We call a 
stationary solution flat if (a,p) = (cr(y),p(y)) and p is a positive constant. Finally, 
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we call the problem ([2]) well-posed if there exist a time T > and a unique solution 
(a,p,p) on [0,T) as defined above. 

The usual way of treating a free boundary problem like © is to transform it 
to a problem on a fixed reference domain. Let f2 := S X (0, 1) with the boundary 
components I\ := § x {i} ~ § for i G {0, 1}. Given any p G C+(§), we introduce 
the map 

p :Ti^TT p , (x',y')^(x',y'p(x')); 

i.e., we will label the coordinates in Q p by (x,y) and the variables in f2 by (x',y') 
with the transformation 

/ y 

x = x and y = y p(x) or x' = x and y' - 



Since p is positive, it follows that 9 P is a C 2 -diffcomorphism. Let u G C(fl p ) and 
v £ C(fl). Then 

0*u := m o 6> p G C(Q) and (0 p )»i; := v o 6T 1 G C(T^). 

We call 6* the pull-back and {Op)* the push-forward operator for the pair (fl, fl p ). 
For v G C 2 (ft) and p G C^(§) wc define 

^(p)« := 9;A[(6 p )*v}, B(p)v := 0;(tr[V(0 p )*v] • i/), 

where tr denotes the trace operator with respect to F p . It is easy to derive the 
explicit formulae 

(3) A(P)V = V X / X > - V X >y> ^ ^ + Vy>y> ~ + ^ ~ ^ ~ Vy'tf — ~ P 



and 



1 . 
Px' + -Vy'\y'=l. 
y' = l P 



The straightforward calculations leading to ([3]) and (j4]) arc omitted for the conve- 
nience of the reader. We conclude that A is uniformly elliptic in f2, as defined in 
[iol ]. Using the fact that the little Holder spaces are Banach algebras, it follows 
that 

A G C 00 (4 +Q (§);£(/i 2+i+Q (fi),/ l l+Q (fi))), i G {0,1}. 

and 

6 G C oc (^ +Q (§);£(/i 2+ ' i+Q (n),/ l 1+i+Q (§))), i G {0,1}. 
In terms of the variables 



r(t):=e; {t) a(t,-) and q(t) := 6* m p(t, ■), 
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the system 



(5) 



is equivalent to 

f A(p(t))r : 
A(p(t))q : 

Pt -■ 

T 
T 

q - 
p - 



-(j,(t - a) 
-B(p(t))q 
02 

1 K 



Po 



in f2, 
in 51, 
on Ti 

on ri 
on r 

on r x 

on r 

for t = 



t > 0, 

t > 0, 

t > 0, 

t > 0, 

t > 0, 

t > 0, 

t > 0, 



0. 



Next, we introduce the following solution operators: The solution of the boundary 
value problem 

A(p)r = t in Q, 
t = cti onr , 
t = (T 2 on Ti 

is denoted as t = lZ(p)(<Ji, 02). For given functions / € /i 1+Q (J7) and fc G h 2+a (E>), 
we write the solution of 



A(p)q 



f hiO, 

on r , 

fc on Ti 



as q = S(p)f + T{p)k. We have 

n{-)(a u a 2 ) e c°°(^ +i+Q (§);^ +l+Q (n)), 

and 



i£{0,l}, 



(5,T) € C°°(^ +Q (§);£(/ l 1+Q (r!),/ l 3+Q (r!)) x £(/i 2+Q (§), /t 2+a (r!))). 

Since the sign of n determines the following results in a very significant way, it 
will be very instructive to remind the reader of the following elementary result: In 
classical differential geometry one considers a local parametrization c(x) = (x, p(x)) 
of the curve T p and has 

= Ms)x^ (s)| = 

Note however that the boundary T p is convex in (x, p(t, x)) with respect to the outer 
normal of f2 p if and only if p(t, x) is concave in x and vice versa. This motivates to 
change the sign of the classical curvature formula so that we will use the identity 

K = -(i + pir 3/2 Pxx 

for the curvature of T p . With this modified sign convention we write 



K r 



V{p)p, 



where 



V(p) = -(l + pl)~ 3/2 d 2 x e C°°(^ +Q (§);/:(^ 4+Q (S),/i 2+Q (§))). 
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3. Well-posedness and the flat stationary solutions 

In this section, we discuss the question of existence and uniqueness of a solution 
to ©. In particular, we are interested in stationary solutions to i.e., solutions 
which do not depend on time. We begin with a proof of the following theorem. In 
fact, the arguments are just a repetition of what is derived in 



121. 



Theorem 1. Given po £ /i^ +Q (§) ; there exists T > such that ^ has a solution 
which is unique in the class C(L> P , T ) x C(D p , T ) x C([0,T) 7 h 3 + a {E>)). Furthermore, 
if Pa 6 h 4 + +a (§), t/ienpe C([0,T),/ l | +Q (§))nC7 1 ([0,T),/ l 1+Q (§)). 

Proof. Using the notation of Section [2J it is easy to see that the transformed sys- 
tem ([5]) can be written in the form 

(6) { Pt + ®{p)P = F(P), t > 0, 
U I P = Po, t = 0, 

where 

$(p) = 7 B(p)T(p)V(p) and F(p)=p,B(p)S(p){TZ(p)(a u a 2 )-a}. 
We find that 

(7) ($,F) G C°°(4 +Q (S);/:(/i 4+Q (S),/i 1+Q (§)) x h 2+a (§)). 

Since 7 > 0, the operator <E>(p) generates, for any p £ /i+ +a (S), a strongly continuous 
analytic semigroup on the space /i 1+Q (§) (cf. Theorem 4.1. in [13 )• An application 
of Amann's local existence, uniqueness and regularity theory for abstract quasilinear 
evolution equations (see, e.g., Theorem 12.1 and Remarks 12.2 in achieves the 
proof. □ 

To find the flat stationary solutions (c*(y),p*{y), p*) of ^ we first solve the 
problem 

f °'l = <T„ 

cr*(0) = ax, 

C7»(/3*) = CT2, 

and obtain the unique solution 

. . . . ._ _ . sinhy 

(8) cr* (y) = t(T2 - a 1 cosh ) -— h a\ cosh y. 

smh p* 

Next, we observe that the unique solution of 

3* = -n(a*-a-), 
= 0, 
= 

is given by 

(9) p*(y) = p— — ai ° 0S P * (y -p^+fj, fa 2 - a*(y) - \a{pl - y 2 ] 

sm p* \ 2 

Since we must demand p'(p*) = 0, we get the condition 

(10) — — -(1 — coshp*) + p* sinhp* = 0. 
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Let a = ^-pL _ suppose 

(11) 0-1,0-2 > o 
to obtain that a > 2. The function 

f a : (0, oo) — > M, i h> q(1 - cosh x) + x sinh x 

clearly satisfies 

lim f a (x) = 0, lim f a (x) = oo, f' a {x) = coshx(a: + (1 — a) tanhx). 

x— >0 ' x— >oo 

By ((TT|) . we conclude that there is a unique p* > with f a (p*) = 0. Now the triple 
(er*,p*,p*) constitutes the unique flat stationary solution of ©. 

The condition (111[) is also reasonable concerning the long-time behavior of our 
model: Let 

Vol(Q p (t)) = / p(i,x)da; 
Jo 

denote the tumor volume. Using the condition p y (x, 0) = and the periodic bound- 
ary conditions for p x , an application of the Gauss-Green Theorem shows that 

Avol(O p (i)) = / pt(t, x) dx = - [ ^ dx = - [ Ap d(x, y). 
at Jo Jdn p (t) ov Jn p (t) 

Applying the maximum principle, we conclude that 

d f 

(12) — Vol(fi p (t)) = /x / (o - a) d{x, y) < n(max{oi,o 2 } - fr)Vol(fi p (*)). 
d * Jn p (t) 

Now (jlip guarantees that the right-hand side of (jT2J) is positive; otherwise we would 
have Vol(f2 p (t)) — > 0, as t — > oo, meaning that the tumor will eventually vanish. 
We have proved the following theorem. 

Theorem 2. Assume that the condition (jlip is satisfied. Then the problem ([2} has 
a unique flat stationary solution (cr*,p*, /)*) which is determined by the formulas 
©, © and (HUD. 



4. The linearization and asymptotic stability 

After studying local solvability and regularity of fully nonlinear equations, the 
second step is to consider the asymptotic behavior and in particular stability of the 
stationary solutions. Recall that a stationary solution u of an autonomous problem 
u'(t) = F(u(t)), t > 0, is called stable if for each e > there is 5 > such that 
for 1 1 i*o — u\\ < 5, we have that \\u(t) — u\\ < e for any t > 0, and the solution 
u = u(t;uo) exists for all i > (denoted as t(uq) = oo). The stationary solution 
u is called asymptotically stable if it is stable and in addition \\u{t) — u\\ — > as 
t — > oo, uniformly for uq in a neighborhood of u. It is said to be unstable if it is 



not stable, [26 1 



In this section, it is our aim to prove that the stationary point (<7»,p*, p*) ob- 
tained in the previous section is asymptotically stable. Therefore, we consider the 
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linearization of ([5]): we plug the ansatz 



= ( P*(y'P*) | + £ I P(t,x',y') 
r{t,x') 

for e > small and with the new unknowns (£, P, r) into ([5]) and compute the 
derivative with respect to e at e = 0. This yields 



(13) 



^x'x' H 


1 v , , 


= 6(cr*)rH 


-s 


in x (0,T), 


Px'x' 


- Xp , , 

1% y v 


= b(p*)r - 


- jiTi 


in x (0,T), 




p y'\ y '=i 


= 




on § x (0,T), 




S|»'=i 


= 




on § x (0,T), 




*V=i 


= ~7?Vx' 




on § x (0,T), 




S|i/'=o 


= 




on S x (0,T), 






= 




on S x (0,T), 




r 


= 7-0, 




f = 0, 



where 



2r 



6( w )r = — v "{y'p*) + »V x /yV(yV,), u e C 2 [0, 1]. 



For a given function r £ h 4+a (E>), we solve the boundary value problem for £ in 
(fTB")) and obtain a unique solution S £ h i+a (£l), which is periodic in x. Substituting 
£ into the second line we get a linear problem for P. Solving the equation for P 
with the corresponding boundary conditions, we obtain a function P £ h 2+a (£l) 
which is also periodic in x; observe that ^r x i x i £ /i 2+Q (S). Let us now introduce an 
operator A £ £(/i 4+Q (§), h 1+a {S)) by setting 



(14) 
Let 



1 BP 

(Ar)(x') = ——(x',l) 
p* oy 



By Q, we have * £ C°°(/i 4+Q (§), h 1+a (S)) and © shows that p t 
solution p of ©. With 



-^(/o) for any 



de 



e=0 



de 



(p* + er) t 



e=0 



and the third equation in (|T3| we see that A = D^(p*). 
As in 12j, we can conclude the following proposition. 

Proposition 3. We have A £ £(ft, 4+a (S), h 1+a (§)), and —A, considered as an 
unbounded operator in h 1+a (E>), generates a strongly continuous analytic semigroup. 



It is our goal to represent the operator A introduced in (|14p as a Fourier multi- 
plication operator. To do so, we use that the functions r £ h i+a (E>), S £ h 4+a (n) 
and P £ h 2+a {VL) are periodic with respect to the variable x' . We thus have the 
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expansions 



r(x') \ / a \ oo / a k b k 
Z{x>,y>) = A (y') \+J2\ Mv') B k {y 
K P(x',y')J \M (y')J k=i\ Mk (y') N^y'), 



( cos(fcx') 
\sin(fca:') 



Proposition 4. The operator A is a Fourier multiplication operator, i.e., given 
r 6 C°°(§) with the Fourier expansion 



r(x') — ao + ^^(afe cos(kx') + b k sm(kx')), 

k=l 

the image Ar has the Fourier expansion 

oo 

{Ar){x') = A a + Afc (a k cos{kx') + b k sm(kx')) 

k=l 

with 

/i(cT 2 cosh p* — ai)\/l + k 2 



A, 



sinh sinh(p*\/l + k 2 ) cosh(/9*fc) 



^cosh(pH,\/l + k 2 ) cosh(/9»fc) - ij 



(15) 



+ jk 2 - flap* - (i 



o~ 2 — o\ cosh p t 



sinh p : 

Moreover, the spectrum cr(A) of A is 

o-(A) = {A fc ; fc e No} 

and o~{A) C K+ for sufficiently large 7. 

Proof. To simplify notation, we will write 
(72 — (7 1 coshp 



fc tanh(yO*fc) + /i((7 — (T2). 



Cl 



j C2 = Ci, C3 = ci coshp, + c 2 sinhp. 



sinh p, 

in the sequel. First, we solve the boundary value problems 



for k <G No, and 



-k 2 A k {y') + ^Al{y') 
A k (0) 
A k (l) 



-k 2 B k (yi) + jrB>i(y') 



B k (0) 
B fc (l) 



A k (y') + a k f k (y'), 

0, 
0. 



B k (y') + b k f k (y'), 

0, 
0. 



for k € N, with 



/fc(2/') = — ( c i sinh(j/'p*) + c 2 cosh^'p*)) - k 2 y'(ci cosh(y'p*) + C2 sinn^'p*))- 
P* 
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A lengthy and somewhat tedious computation shows that 



A k {y ) = flfeCi y cosh(y p„) — coshp* 

\ smh(p„Vl + «) / 

sinh^yWl + fc 2 ) . \ 

+a fc c 2 y smh(y p„) — smhp* , 

\ smh(p*v 1 + fc ) J 

and B k (y') is obtained from Afc(y') by exchanging a k with fe^. The next task is to 
solve the boundary value problems 

-k 2 M k (y') + jl M 'k(v') = -vA k (yi) + a k g k (yi), 
* M k (0) = 0, 
Mfc(l) - 7 fc 2 a fc , 



for fc g No, and 



-fc 2 A,(j/) + £iW) = -l*BkW) + b k g k (3/), 
* ^(0) - 0, 



JV fc (l) = 7 fc 2 & fe , 



for fc e N, with 



2 

5fe(y') = — K° - a sinh(y'p*) ~ c 2 cosh(y'p*)) 
P* 

-k 2 p(ciy' + crp*{y') 2 - c 2 y sinh(?/p*) - Ciy' cosh(j/'p*)). 
Again, it is straightforward to derive the solutions 

M (y') = - 3 (y'p* - p* - sinb/y'p*)) - aofrirp* - a pci 
smh p* 

+a pciy' + a pap*(y') 2 - a p(ciy' cosh(y'p*) + c 2j/' sinh(y'p*)) 

and, for fc 7^ 0, 



M fe (j/) = - - ^^^^ fy ( sinh(y'M) - tanh(p.fc) cosh^M)) 

sinh(y'p»\/l + fc 2 ) , ,2 - . cosh(y'p*&) 

+a fc pc 3 — — +a fc (7fc - pcrp* - pxi) — — — 

smh(p*vl + fc ) cosh(p*fc) 

+a k pc 1 y' + a k p,ap*{y') 2 - ak^iy' cosh(y'p*) + c 2 y' sinh(j/'p*))i 

the A& are obtained from by replacing a k with 6^,. By the definition of A, we 
see that A is a Fourier multiplication operator with 

\ k = — M£(l). 
afep* 

Since /i 4+Q (§) is compactly embedded into /i 1+t *(§) and the resolvent set p(A) is 
non-empty, the resolvent operator (A — A) -1 is compact for any A ^ &(A). Hence 
the spectrum a(A) of A consists entirely of eigenvalues. Our explicit computations 
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show that, for all k £ No, 



pc 3 Vl + k 2 nc 3 vl 

Afc = 



sinh(p*Vl + k 2 ) cosh(p*fc) tanh(p»\/l + k 2 ) 
+('jk 2 — \iap* — ^ci)fctanh(p*fc) 

+2pa + — (ci — C3) — fx(ci sinhp* + C2 coshp*), 
P* 

and, by our definitions, Afc is as specified in (|f 51) . Concerning positivity, we first 
concentrate our attention on An.. Using once again Eq. (|10p . we obtain 



A = p,a 1 - 



sinh p» 

and written in this form, An, is clearly positive. For k ^ 0, fix some 70 > and let 
Afc (70) be as in (fT5|) . Since 

Afc (70) n c , 

Ptanhfo k) = 70 + nfc ~* ° ^ °°' 

we see that there exists fco € N such that Afc (70) > for all k > fco- Since 
Afc (7) > Afc (70) for 7 > 70, it follows that Afc (7) > for all k > k and 7 > 70- It 
is obvious that we can choose 7 > 70 so large that also Ai(7), . . . , Afc (7) > 0. Let 
A 7^ be the collection of all 7 > such that Afc (7) > for all k € N. We have 
shown that a (A) C R.+ for any 76 A and this achieves the proof. □ 
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Figure 2: The spectrum of the operator A, for \i = a = 7 = 1, u\ = 2, 01 = 3. 



In fact, not the operator A but —A will determine the stability properties of the 
problem @. We will employ the following stability theorem. 



Theorem 5 (see |26(). Let X be a Banach space and let T: D(T) C X — > X be a 
linear sectorial operator such that the graph norm of T is equivalent to the norm of 
D. Assume furthermore that 

sup{Rc A, A G er(T)} = -u < 0. 
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Let O be a neighborhood of zero in D and let G : O — > X be a C 1 function with 
locally Lipschitz continuous derivative such that G(0) = and G'(0) = 0. Fix 
Li € [0, wo)- Then there exist r > and M > such that for each uq € -6(0, r) C D 
the solution u(t) to the problem 

u'(t) = Tu(t) + G(u(t)), t > 0, u(0) = u 0) 

satisfies t(uq) = oo and 

\\u(t)\\ D + \\u'(t)\\ x < Me-^\\u \\D, t > 0. 

As for the model in [lij , we show that the steady state of 

(16) p* + *(p) = 0, p(0)=p o , 

is asymptotically stable under small perturbations belonging to /i + +Q (§). By defi- 
nition, = $(p)p - F(p), for p e fr+ +a (§). Setting G(r) := + p*) - Ar, we 
have that 

G g C 00 {h 4 + a (§),h 1+a (S)) 

and observe that 

G(0) = #(p*) = 0, L>G(0) = - A = 0. 

Theorem [5] implies that the solution r to 

r t = -Ar - G(r) 

is asymptotically stable: There are constants hi, e, K such that if ro € /i+ +Q (§) 
satisfies ||t'o||c 4 +° < £ then the solution r exists globally and satisfies 

||r(f)||(7<+- < Acxp(-^)||r ||c 4 +°, t > 0. 

Letting r(t) = p(t) — p» for t > 0, this in turn shows that if ||po — P*!lc 4 +° < £ i 
then the solution to (|16p exists globally and satisfies 

(17) \\p(t) - p*\\c*+o < K exp(-uit)\\p -p*\\ c *+o, t>0. 

Next let us consider the dynamical behavior of a and p. By definition, we have 

ct* = 7£(p*)(o-i,ct 2 ) an d c(t) = 72.(p(t))(ai, ct 2 ). 
By the mean value theorem, there exists a constant G such that 

IKf)-<7*|| C 4 +Q = ||[ft(p(f))-ft(pO]fa.fra)||c*+<> < c\\ P (t) -p4c*+° 

for any t > 0. Combining this with estimate (|17[) . we get 

|| cr(i) — cr*||c4+c« < ATexp(— wi). 

A corresponding estimate for p can be obtained similarly. We have proved the 
following theorem. 

Theorem 6. Let a±,&2 > o > 0, p > and 7 g A oe given. Then the flat 
stationary solution defined by (J8J), ([9]) and ()10[) is asymptotically stable: There are 
positive constants h>, K and e such that if ||po ~~ P*||c 4 + Q < e 

IK*) - o-*llc 4 +° + UK*) -p*\\c*+<> + \\p(t) - P*llc 4 +° < A'cxp(-cjt), 
/or an?/ t > 0. 
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